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We report a measurement of the mass of the top quark in lepton+jets final states of pp —¥ti 
data corresponding to 2.6 fb _1 of integrated luminosity collected by the DO experiment at the 
Fermilab Tevatron Collider. A matrix-element method is developed that combines an in situ jet 
energy calibration with our standard jet energy scale derived from studies of 7+jet and dijet events. 
We then implement a flavor-dependent jet response correction through a novel approach. This 
method is used to measure a top-quark mass of m t = 176.01 ± 1.64 GeV. Combining this result 
with our previous result obtained on an independent data set, we measure a top-quark mass of 
mt = 174.94 ± 1.49 GeV for a total integrated luminosity of 3.6 fb _ . 

PACS numbers: 14.65.Ha 



I. INTRODUCTION 

The observation of the top quark in 1995 P, Q con- 
firmed the existence of the six quarks in three genera- 
tions of fermions expected in the standard model (SM) 
of particle interactions. Because of its mass, the lifetime 
of the top quark is much shorter than the time-scale of 
hadronization. The top quark can therefore decay before 
interacting, making it the only quark whose characteris- 
tics can be studied in isolation. The large mass of the top 
quark (mt), corresponding to a Yukawa coupling to the 
Higgs boson equal to 1 within the current uncertainties, 
suggests a special role for the top quark in the breaking of 
electroweak symmetry. It is therefore not surprising that 
the precise determination of the mass of the top quark has 
received great attention. The interest in the top-quark 
mass also arises from the constraint imposed on the mass 



of the Higgs boson, m# , from the relationship among the 
values of m t , run , and the SM radiative corrections to the 
mass of the W boson Q. A precise measurement of m t 
also provides a useful constraint on contributions from 
physics beyond the standard model 

The statistical uncertainty on the world average value 
of m t is 0.3%, and the accuracy of the measurement of 
m t is now dominated by systematic uncertainties . The 
main systematic contributions arise from uncertainties on 
the jet energy calibration and on the Monte Carlo (MC) 
simulation of ti events. 

We present a new measurement of the mass of the top 
quark based on 2.6 fb _1 of integrated luminosity from 
pp collisions at ^fs = 1.96 TeV, collected with the DO 
detector at the Fermilab Tevatron Collider. The analy- 
sis focuses on ti events identified in lepton+jets (£+jets) 
final states (with £ representing either an electron or a 
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muon) Q, in which the top and antitop quark are as- 
sumed to decay into a W boson and b quark Q, with 
one of the W bosons in the W + W~bb final system de- 
caying via W — > tv into a lepton and neutrino and the 
other via W — > qq' into two quarks, and all four quarks 
(qq'bb) hadronizing into jets. Such events are character- 
ized by an isolated electron or muon with large transverse 
momentum (pt), an undetected neutrino that causes a 
large imbalance in transverse momentum, and four high- 
Pt jets. In selecting candidate events, we exploit this 
distinct signature, which helps distinguish these events 
from background. 

Compared to the previous measurement based on data 
corresponding to 1 fb _1 of integrated luminosity Q, we 
use a larger data set and an improved evaluation of sys- 
tematic uncertainties. The analysis uses the same matrix 
element (ME) analysis technique, with an in situ jet en- 
ergy calibration based on constraining the invariant mass 
of the two jets from the decay of the W boson to the world 
average value of M w = 80.4 GeV @|. As in the previ- 
ous measurement, the standard jet energy scale (JES), 
derived from 7+jet and dijet data samples, is used as an 
additional constraint, and implemented through a Gaus- 
sian prior on its absolute value and uncertainty. A major 
improvement in this new measurement is the significant 
reduction of the uncertainty associated with the model- 
ing of differences in the calorimeter response to 6-quark 
and light-quark jets originating from the introduction of 
a new flavor-dependent jet energy response correction. 

This measurement, like all direct measurements of rrit, 
relies on MC ti events for absolute calibration. It is 
therefore important to understand the precise definition 
of the input mass mf en in MC ti event generators, such 
as ALPGEN [13 and pythia used to calibrate the 
direct measurements. Although mf cn is not well defined 
in leading order (LO) generators that use parton show- 
ers to model higher-order effects and hadronization, it 
has been argued that mf en should be viewed as being 
close to the pole mass [12]. In Ref. [l3[, the DO Col- 
laboration has extracted mt from a comparison of the 
measured ti production cross section with predictions 
from higher-order quantum chromodynamics (QCD), by 
equating mf cn both with the pole mass (mj? ole ) and with 
the MS mass (m^). The extracted m t , under the as- 
sumption mf en = m^ olc , is found to agree with the aver- 
age value of mt from the Tevatron, while the mt extracted 
assuming mf on = m^ s is found to be different from the 
average value of mt. These results favor the pole mass 
interpretation of mf cn . 

This paper is arranged as follows. A brief description 
of the DO detector is given in Sec. HU which is followed 
by a discussion of the selection and reconstruction of the 
physical objects in this analysis in Sec. IIIII Section II VI 
summarizes the MC samples used to simulate the events 
of interest, and Sec. [V] discusses the technique used to 
extract the value of m t . This is followed by a description 
of the calibration of the response of the analysis method 



in Sec. I VII and a discussion of the flavor-dependent jet 
response correction used to bring the simulation of the 
calorimeter response to jets into agreement with data in 
Sec. I VIII The result of the calibration is applied to the 
data in Sec. IVIII1 where the measured value of m* and 
its statistical uncertainty are also presented. Section IIXI 
describes the evaluation of systematic uncertainties and 
the final result is given in Sec. |Xj We combine this new 
measurement in Sec. IXII with an updated version of that 
from Ref. Q in which the flavor-dependent jet response 
correction mentioned above has been applied and the sys- 
tematic uncertainties have been updated. 



II. THE DO DETECTOR 

The DO detector consists primarily of a magnetic 
central tracking system, calorimetry, and a muon sys- 
tem. The central tracking system comprises a silicon mi- 
crostrip tracker (SMT) and a central fiber tracker (CFT), 
both located within a 1.9 T superconducting solenoidal 
magnet [H. The SMT [H has w 800,000 individual 
strips, with typical pitch of 50 — 80 /im, and a design op- 
timized for track and vertex finding at |ry| < 2.5, where 
the pseudorapidity rj = — In [tan(#/2)], and 9 is the po- 
lar angle with respect to the proton beam direction rel- 
ative to the center of the detector. The system has a 
six-barrel longitudinal structure, each with a set of four 
layers arranged axially around the beam pipe, and in- 
terspersed with 16 radial disks. In 2006, a fifth layer, 
referred to as Layer 0, was installed close to the beam 
pipe [lj| [l7|- The CFT has eight thin coaxial barrels, 
each supporting two doublets of overlapping scintillating 
fibers of 0.835 mm diameter, one doublet being parallel 
to the collision axis, and the other alternating by ±3° rel- 
ative to the axis. Light signals are transferred via clear 
fibers to solid-state photon counters (VLPCs) that have 
~ 80% quantum efficiency. 

Central and forward preshower detectors, located just 
outside of the superconducting coil (in front of the 
calorimetry) , are constructed of several layers of extruded 
triangular scintillator strips that are read out using 
wavelength-shifting fibers and the VLPC. These detec- 
tors provide initial sampling of electromagnetic showers, 
and thereby help distinguish incident photons from elec- 
trons. The next layer of detection involves three liquid- 
argon/uranium calorimeters: a central section (CC) cov- 
ering \rj\ up to ss 1.1, and two end calorimeters (EC) 
that extend coverage to \r]\ ss 4.2, all housed in sepa- 
rate cryostats. The electromagnetic (EM) section of the 
calorimeter is segmented into four layers, with transverse 
segmentation of the cells in pseuodorapidity and azimuth 
of A77 x Acf) = 0.1 x 0.1, except for the third layer, where 
the segmentation is 0.05 x 0.05. The hadronic portion 
of the calorimeter is located after the EM sections, and 
consists of fine hadron-sampling layers, followed by more 
coarse hadronic layers. In addition, scintillators between 
the CC and EC cryostats provide sampling of developing 
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showers for 1.1 < |nj < 1-4. 

A muon system [18| is located beyond the calorimetry, 
and consists of a layer of tracking detectors and scintil- 
lation trigger counters before 1.9 T toroids, followed by 
two similar layers after the toroids. Tracking for < 1 
relies on 10 cm wide drift tubes, while 1 cm mini-drift 
tubes are used for 1 < \r]\ < 2. 

Luminosity is measured using plastic scintillator ar- 
rays located in front of the EC cryostats, covering 2.7 < 
|?7| < 4.4. The trigger and data acquisition systems are 
designed to accommodate the high instantaneous lumi- 
nosities of the Tevatron [TJ, Based on preliminary 
information from tracking, calorimetry, and muon sys- 
tems, the output of the first level of the trigger is used to 
limit the rate for accepted events to rj 2 kHz. At the next 
trigger stage, with more refined information, the rate is 
reduced further to rj 1 kHz. These first two levels of trig- 
gering rely mainly on hardware and firmware. The third 
and final level of the trigger, with access to all of the 
event information, uses software algorithms and a com- 
puting farm, and reduces the output rate to w 100 Hz, 
which is written to tape. 



III. OBJECT RECONSTRUCTION AND EVENT 
SELECTION 

In the following sections, we summarize how the physi- 
cal objects in data and MC events are reconstructed from 
information in the detector and the criteria applied to 
these objects to select the £+jets ti candidate events. 

A. Object Reconstruction 

This section describes the reconstruction of electrons, 
muons, missing transverse momentum, and jets, and the 
identification of b jets. 

1. Identification of Electrons 

Electron candidates are defined by narrow clusters 
of energy deposited in towers of the electromagnetic 
calorimeter located within a cone of radius TZ = 
J (A77) 2 + (A(p) 2 = 0.2, where Arj is the pseudorapid- 
ity, and A</> is the azimuthal angle of each cluster rel- 
ative to the seed cluster. At least 90% of the total 
energy measured within this cone is required to be lo- 
cated within the electromagnetic section to be consis- 
tent with expectations for electromagnetic showers. Iso- 
lation from energy deposited by hadrons is imposed by 
requiring (E to t - Eem)/Eem < 0.15, where E to t (Eem) 
is the total (electromagnetic) energy in a cone of radius 
TZ = 0.4 (TZ = 0.2). Candidate electrons are required 
to have longitudinal and transverse shower profiles com- 
patible with those of electromagnetic showers and to be 
spatially matched to a track reconstructed in the central 



tracking system. Electron candidates meeting these cri- 
teria are referred to as loose electrons. Finally, (i) the 
value of a multivariable likelihood discriminant based on 
tracking system and calorimeter information is required 
to be consistent with that for an electron, and (ii) a neu- 
ral network, trained using information from the tracking 
system, calorimeter, and central preshower detector is 
used to further reject background from jets misidentified 
as electrons. Electron candidates meeting these criteria 
are referred to as tight electrons, and are those used to 
obtain the final selection. 



2. Identification of Muons 

Muons are identified by requiring a minimum number 
of wire and scintillator hits on both sides of the toroidal 
magnets in the muon detector [18| . Cosmic ray back- 
ground is rejected by requiring scintillator signals consis- 
tent in time with muons originating from the pp collision. 
Tracks in the muon system are required to match a re- 
constructed track in the central tracker having a small 
impact parameter with respect to the pp interaction ver- 
tex (PV) to reject muons from cosmic rays and decays in 
flight of kaons and pions. Muon candidates must also 
be isolated from jets with p^ > 15 GeV by requring 
a separation in r\ — (j) space between the muon and jet 
of A7?.(/x,jet) > 0.5 [20| . Candidates satisfying these 
requirements are referred to as loosely isolated muons. 
The following two variables are used to impose additional 
isolation requirements: E^^ d is defined as the ratio of 
calorimeter energy within an annulus of 0.1 < TZ < 0.4 
around the muon direction to the pr of the muon; cone 
is defined as the ratio of the total px of all tracks within 
a cone of TZ = 0.5, excluding the muon, to the pr of the 
muon. Muon candidates meeting all the requirements 
above that satisfy E^ l o cd < 0.12 and p^cone < °- 12 are 
referred to as veto muons. Further tightening these re- 
quirements to ££ c a t cd < 0.08 and p s ^* e < 0.06 selects 
candidates referred to as tightly isolated muons. 

3. Measurement of the Imbalance in Transverse Momentum 

We use the conservation of momentum to measure the 
momentum imbalance in the transverse plane (^t)- From 
that, we infer the presence of the neutrino. The j>T is 
determined from the vector sum of the energies of all 
cells in the electromagnetic and hadronic calorimeters. 
Subsequent energy corrections applied to reconstructed 
objects such as jets and muons are also propagated to 
the missing transverse momentum. 

4- Identification of Jets 

Jet candidates are reconstructed using the iterative 
midpoint cone algorithm with a cone radius of TZ = 
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0.5 [2l|. Only calorimeter cells with energies that are 
2.5 standard deviations above the mean of the noise dis- 
tribution are considered in the reconstruction. Isolated 
cells with energies less than 4 standard deviations above 
the mean of the noise distribution are also discarded. 
Among the jet candidates with pr > 8 GeV, the follow- 
ing selection criteria are imposed. The electromagnetic 
fraction of the jet energy is required to be below 0.95 
to reject electrons and above 0.05 to suppress jets domi- 
nated by noise from the hadronic part of the calorimeter. 
Jets with a large fraction of their energy deposited in the 
coarse hadronic layers of the calorimeter are rejected to 
suppress jets dominated by noise typical for those layers. 
To minimize background from jet candidates arising from 
noise in the precision readout of the calorimeter, confir- 
mation from the readout system of the first level trigger 
is required for reconstructed jets. Jets matched to loose 
electrons with p T > 20 GeV and Aft(ejet) < 0.5 [U 
are also rejected. Energies of jets containing muons are 
corrected with the measured muon momentum after ac- 
counting for the typical energy deposited by a minimum 
ionizing particle. 

The energy of a reconstructed jet is corrected, on aver- 
age, to that of a particle jet [22| containing the final-state 
particles within a cone of radius 1Z = 0.5 corresponding 
to the reconstructed jet. The first step involves the sub- 
traction of the offset energy due to calorimeter noise and 
contributions from previous and following beam crossings 
and multiple interactions within the same beam crossing. 
This is followed by an absolute response correction de- 
termined from 7+jet events and a relative ^-dependent 
correction based on 7+jet and dijet events. Finally, a 
showering correction is applied to account for the lateral 
leakage of energy across the jet cone boundary. 



5. Identification of b jets 

The lifetime of the b quark, unlike that of the top 
quark, is far longer than the time-scale for hadronization. 
This means that, during QCD evolution, the b quark can 
form short lived b hadrons that travel > 1 mm before de- 
caying through the weak interaction. We identify the b 
jets among the candidates satisfying the jet selection cri- 
teria described in the previous section by using a neural 
network (NN) 6-tagging algorithm that selects jets with 
displaced vertices and tracks relative to the P V |23| ■ The 
NN tagger is based on nine input variables that can be 
separated into two categories. The first category is re- 
lated to the reconstructed secondary vertex and includes 
the vertex quality, the number of associated tracks, the 
invariant mass of the vertex, the number of secondary 
vertices reconstructed within the jet, the spatial separa- 
tion between the jet axis and the position vector of the 
secondary vertex relative to the PV, and the length of 
the flight path projected on the transverse plane divided 
by its uncertainty (which provides a measure of the decay 
length significance in terms of standard deviations) . The 



second category relies only on the characteristics of the 
tracks within the jets such as impact parameters, trans- 
verse momentum, and track quality. The &-jet candidates 
are also required to have at least two good quality tracks 
originating from the PV. The tagging efficiency for b jets 
is m 65% for a misidentification rate of ss 3% for u, d, s 
quark, or gluon jets (24l |. 

B. Event Selection 

The data sample used in this analysis was collected 
with the DO detector at the Tevatron between June 2006 
and June 2008, and corresponds to an integrated lu- 
minosity of 2.6 fb _1 . The selected events must satisfy 
a single-lepton trigger, requiring a high px electron or 
muon, or a lepton+jets trigger, requiring a lower-pT elec- 
tron or muon accompanied by a jet. Events are required 
to have at least one PV with > 2 tracks reconstructed 
within the fiducial region of the SMT. We require ex- 
actly four jets with |?7| < 2.5, with the leading (highest 
Pt) jet having pt > 40 GeV, and the other jets pr > 20 
GeV. Leptons are required to originate from within 1 cm 
of the PV in the coordinate along the beam line. Ex- 
actly one tight electron (or tightly isolated muon) with 
Pt > 20 GeV and \r)\ < 1.1 (|?7| < 2) is also required. 
Electron+jets events containing a second tight electron 
with pt > 15 GeV and |t;| < 2.5 or a veto muon with 
Pt > 15 GeV and \r)\ < 2 are rejected. Muon+jets 
events containing a second muon that is a veto muon with 
Pt > 15 GeV and \r)\ < 2 or a tight electron with pt > 15 
GeV and \r]\ < 2.5 are rejected. The missing transverse 
momentum is required to satisfy ■pr > 20 GeV (J>t > 25 
GeV) for e+jets (/x+jets) events. Multijet background, 
typically arising from mismeasurement of lepton or jet 
energies, is suppressed by requiring a minimal azimuthal 
separation between the lepton direction and the j>T vec- 
tor with A(j)(e,pT) > 0.7-7T — 0.045 • j>T for electrons and 
A<fi(fj,,pT) > 2.1 — 0.035 • j>T for muons, with j>T in GeV 
and A4>(£,Pt) = \4>i ~ 4>itr\- Any ^+jets events with an 
invariant masss, m MM , of the isolated muon and a second 
muon (with pt > 15 GeV and even lower quality require- 
ments than a loosely isolated muon) of 70 < < 110 
GeV are rejected in order to suppress Z(—¥ /i/i)+jets 
events. The data sample satisfying the above criteria 
consists of 825 e+jets and 737 /x+jets events. We fur- 
ther require at least one jet to be identified as a b jet, 
which yields the final data samples of 312 e+jets and 303 
/i+jets events. 



IV. MONTE CARLO SAMPLES 

The MC events used to model the it signal and the 
W+jets background needed for the calibration of the 
measurement (described in Sec. IVI|) are generated using 
ALPGEN [Toll to simulate the hard-scattering process and 
pythia 11 ll to simulate hadronization and shower evo- 
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lution. The MLM matching scheme (25J is employed to 
avoid overlaps between components of the event belong- 
ing to the hard process, implemented through a matrix 
element, and parton evolution (showering) into jets. The 
W+jets background samples are divided into two cate- 
gories: (i) W + lp and (ii) W + (cc, Ip) and W + (66, Ip), 
where Ip (light partons) denotes u, d, ,s-quarks, or glu- 
ons. Although the individual processes are produced with 
ALPGEN which is a LO generator, the relative contri- 
butions between the two categories are determined us- 
ing next-to-leading order (NLO) calculations, with next- 
to-leading logarithmic (NLL) corrections based on the 
mcfm MC generator (2g. The MC samples used to de- 
rive jet transfer functions that correlate jet energies with 
those of partons in the it events (described in Sec. IV B 1J 
are generated using pythia to simulate both the hard- 
scattering process and the subsequent hadronization and 
shower evolution for the events. All MC samples are 
generated with CTEQ6L1 parton distribution functions 
(PDFs) and passed through a full GEANT3-based [H] 
simulation of the DO detector. To simulate the effects 
from additional pp interactions, events with no trigger 
requirements selected from random pp crossings in the 
collider data having the same instantaneous luminosity 
profile as the data are overlayed on the fully simulated 
MC events. This is then followed by the same reconstruc- 
tion and analysis chain as applied to data. 



V. METHOD OF ANALYSIS 

The top-quark mass is measured using all kinematic 
information with a likelihood technique based on prob- 
ability densities (PDs) constructed, for each individual 
event, from matrix elements of the processes contribut- 
ing to the observed final state. This analysis technique, 
referred to as the matrix element (ME) method (2^, is 
described below. 



A. Matrix Element Method 

If the processes contributing to an observed event do 
not interfere, the total PD for observing a given event is 
the sum of all contributing probabilities for that specific 
final state. Assuming that ti and W+jets production are 
the only two contributions, the PD for observing each 
event is given in terms of the top-quark mass rat, the 
jet energy scale factor fcjES dividing the energies of all 
jets, and the fractions of it signal (/) and of VF+jets 



background (1 — /) in the data by: 

Povt = A(x)[fP sig (x; m t , fcjEs) 

+(i-/)iW(s;fcjEs)]. (!) 

where x represents the measured jet and lepton energies 
and angles; A(x), which depends only on x, accounts for 
the geometrical acceptance and efficiencies; and P s i g and 
Pbkg are the PDs for ti and VF+jets production, respec- 
tively. For events satisfying Pbkg 3> Psig, the relative 
contribution of P slg to P ovt is negligible and has minimal 
influence on the determination of mt- Multijet events 
satisfy this condition and can therefore be represented 
by Pbkg, as the event kinematics are far closer to VF+jets 
than to ti production. 

Due to the finite detector resolution and the hadroniza- 
tion process, the measured set x for the observed events 
will not, in general, be identical to the corresponding set 
y of the original final-state partons and the relationship 
between x and y is described by a transfer function. In 
addition, the initial partons carry momenta q\ and qi in 
the colliding p and p. To account for this complication, 
P s ig and Pbkg must be integrated over all parton states 
contributing to the observed set x. This involves a convo- 
lution of the partonic differential cross section da(y) with 
the PDFs and a transfer function W(x,y; &;jes) that re- 
lates x and y: 

P(x,a) = — [ V do(y,a)dqidq 2 f(qi)f(q2) 

flavors 

x W(x,y; fcjEs), (2) 

where a represents the parameters to be determined in 
the analysis, the sum runs over all possible initial-state 
parton flavors, and f(qi) are CTEQ6L1 PDFs for finding 
a parton of a given flavor and longitudinal momentum 
fraction qi in the p or p. Detector resolution is taken 
into account in W(x,y; £;jes), representing the probabil- 
ity density for the measured set x to have arisen from the 
partonic set y. Dividing by the total observed cross sec- 
tion for the process, a(a), ensures P(x; a) is normalized 
to unity. 

The differential cross section term for P s i g is calcu- 
lated using the LO ME of the quark-antiquark annihi- 
lation process (Ai t j). A total of 24 integration vari- 
ables are associated with the two initial state partons and 
the six particles in the final state. Since the angles for 
the four jets and the charged lepton are sufficiently well 
measured, the angular resolution terms in W(x,y; /cjes) 
can be approximated by Dirac 5-functions. Integrating 
over these and four more (5-functions that impose energy- 
momentum conservation leaves 10 integrals to evaluate 
the probability density that represents the it production 
probability for a given mt and /cjes (30| : 
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P sig = -^X> / d P dm\dMldm\dMld Pl dqfdq\d<%dql £ |7W^| 2 L (gl)/ fe) $ 6 W(s,y; fcjE S ),(3) 

flavors,!/ V^a^ltff) 2 -™^™^ 
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(muon), and the energies of the four partons producing 
the jets, summing over the 24 jet permutations and all 
neutrino solutions. 

The top-quark mass is extracted from n events with a 
measured set of variables x = (x±, X2, x n ) through a 
likelihood function for individual event probabilities P evt 
according to 

n 

L(x; m t , fcjEs, /) = \\ P C vt(xf, m t , fcj E s, /)• (4) 



Figure 1. (color online) The mean acceptance as a function 
of mt and /cjes for the (a) e+jets and (b) /x+jets channels. 

where, in addition to the CTEQ6L1 PDF given by 
the / (qi) also include parameterizations of the probabil- 
ity distributions for the transverse momenta qf ,y of the 
initial-state partons in PYTHIA The masses of the 

initial-state partons are given by m qi , and $6 includes 
the 6-body phase-space factor and other constants. The 
first sum is over all 24 jet permutations, each carrying a 
weight Wi, which is the product of four jet weights. The 
weight for a 6-tagged jet with a given pt and rj is the av- 
erage tagging efficiency e a {pT,vi) for a given parton hy- 
pothesis a (= b, c, light q, or gluon). The weight for a jet 
that is not 6-tagged is 1 — e a (pT, r/). The second sum in- 
cludes up to eight solutions for neutrino kinematics, and 
conservation of transverse momentum used to calculate 
the transverse momentum of the neutrino. The parame- 
ter p represents the fraction of the energy carried by one 
of the quarks from the W — > qq' decay. The masses of 
the two W bosons (M\,M2) and of the pair of top quarks 
(mi, 7712) are chosen as integration variables because of 
computational efficiency related to the four Breit-Wigner 
mass terms that make the ME negligible everywhere 
except at the mass peaks. The energy (the curvature 
1/pr) of the electron (muon) is defined by pt. The in- 
tegration over qi involves only transverse components. 
W{x, y; fcjEs) is the product of five terms for the four jets 
and one charged lepton, described below. The normal- 
ization a" bs = / A(x)P s i g dx = a u (m t ) x (A{m t , /cjes)) 
is calculated from the product of the total cross section 
corresponding to the ME used and the mean acceptance 
for events whose dependencies on mt and /cjes are deter- 
mined from MC events. The mean acceptance is shown 
in Fig. [1] as a function of m t for different values of fcjEs- 

The differential cross section in Pbkg is calculated using 
the W+4: jets matrix elements from the VECBOS HH MC 
program. The initial-state partons are assumed to have 
no transverse momenta. The integration is performed 
over the W boson mass, the energy (1/pr) of the electron 



For every assumed pair of (m tl fcjss) values, the value 
of / bcst that maximizes the likelihood is determined. 
To obtain the best estimate of m t and fcjES; the two- 
dimensional likelihood: 

L (x; mt, fcjEs) = L [x; m t , fc JE s, f hcst (m t , fcj E s)] (5) 

is projected onto the mt and fcjss axes according to 

L (x; m t ) = J L(x; mt, fc JE s) G (fc JE s) dfcjES (6) 

and 

L (x; fcj E s) = / L (x; m t , fcjEs) dmt, (7) 

using Simpson's rule [32l |. where the prior probability 
distribution G(fcjEs) is a Gaussian function centered 
at fcjEs = 1 with standard deviation (sd) 0.02 deter- 
mined from the mean of the fractional uncertainty of 
the standard jet energy scale corrections applied to all 
jets in the MC samples used in this analysis. The 
best estimates and the uncertainties on the mass of 
the top quark and the jet energy scale are then ex- 
tracted using the mean and the RMS of L (x; m t ) and 
L (x; fcjEs); respectively. The mean is calculated from 
a = J aL(x; a) da/ J L(i; a)da and the RMS from 
cr 2 (a) = / (a — a) 2 L(x; a) da/ J L(x; a) da, where a 
corresponds to m t or /cjes, also using Simpson's rule. 

The fit parameter ZcjeSi associated with the in situ jet 
energy calibration, has the effect of rescaling the energies 
of all the jets, and thereby the 2-jet invariant mass of the 
hadronically decaying W boson, with the jet energy scale 
factor /cjes ■ The presence of the Breit-Wigner mass term 
associated with the hadronically decaying W boson in the 
ME of Eq. ([3]) maximizes the likelihood in Eq. |4|) when 
the 2-jet invariant mass coincides with the Breit-Wigner 
pole fixed at the world average of M w = 80.4 GeV @. 
The additional constraint to the standard scale derived 
from 7+jet and dijet samples is applied through the prior 
probability distribution G(/cjes) in Eq. 
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Table I. Transfer function parameters for light quarks (cn in 
GeV). 
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Detector Resolution 





In this section, we describe the parameterizations for 
the jet and electron energy and muon px resolutions used 
in the transfer function W(x,y; fcjEs) which is the prod- 
uct of four jet transfer functions for a given jet permuta- 
tion and an electron or muon transfer function. 



1. Parameterization of Jet Energy Resolution 

The transfer function for jets, Wj e t(E x , E y ; fcjEs), rep- 
resents the probability that a measured jet energy E x in 
the detector corresponds to a parent quark of energy E y . 
It is parameterized in terms of a double Gaussian func- 
tion whose means and widths are dependent on E y . For 
the case /cjes = 1, it is given by 



U>, (/.', . /.'.,: ^es = 1) 



1 



/ 2tt(P2 +P3P5) 

^? 



P3e 



-Eb)-P4] 




,(8) 



where the pi are functions of the quark energy for quark 
i and are parameterized as linear functions of the E y : 
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Figure 2. (color online) Transfer functions for fcjES =1 light- 
quark jets as a function of measured jet energy for different 
parton energies in n regions: (a) \n\ < 0.5, (b) 0.5 < \r)\ < 1.0, 
(c) 1.0 < \n\ < 1.5, and (d) 1.5 < \r)\ < 2.5. 



for each event with respect to a; and bi. A different set 
of parameters is derived (i) for three varieties of quarks: 
light quarks (u, d, s, c), b quarks with a soft muon tag 
in the jet |33j, and all other b quarks, and (ii) for four 
7? regions: \r)\ < 0.5, 0.5 < |ry| < 1.0, 1.0 < |?7| < 1.5, 
and 1.5 < < 2.5, to minimize possible effects due to 
non-uniform calorimeter response. The values for these 
parameters are shown in Tables U and |TT] for light-quark 
and 6-quark jets, respectively. Figure [2] illustrates the 
transfer functions for light-quark jets as a function of E x 
for different values of E y . In Fig. [31 we compare the 2- 
jet and 3-jet invariant mass distributions for two types 
of pythia ti £+jets events: (i) parton level events with 
jet energies smeared using the transfer functions and (ii) 
fully simulated events where all four reconstructed jets 
are matched to partons with A7?.(parton, jet) < 0.5. The 
2-jet (3-jet) invariant masses are calculated using the two 
light-quark jets (all three jets) from the hadronic branch 
of the ti £+jets events and correspond to the W boson 
(top-quark) mass. The overlaid distributions in Fig. [3] 
indicate that the jet transfer functions describe the jet 
resolutions well. 

For A;jes 7^ lj the jet transfer function is changed to 



The parameters Oj and bi are determined from fully sim- 
ulated ti events, following all jet energy corrections and 
smearing to match resolutions in data. These events are 
generated with pythia at nine values of the top-quark 
mass ranging from 155 to 195 GeV in 5 GeV intervals. 
The parton and jet energies are used in an unbinned like- 
lihood fit that minimizes the product of the Wj e t terms 



Wj e t (E x ,E y ; fcjEs) = 



Hi 



( E * 



,E y ; 1 



«JES 



(10) 



where the /cjes factor in the denominator preserves the 
normalization J Wj ct {E x , E y ; kj ES )dE x = 1. 
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Table II. Transfer function parameters for 6-quark jets without and with a muon within the jet cone (a; in GeV). 
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6-quark jets with a muon within the jet cone 
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Figure 3. (color online) Comparison of (a) 2-jet and (b) 3- 
jet invariant mass distributions for parton-level ti MC events 
with energies smeared using the transfer functions (open his- 
togram) and fully simulated tt MC events with all four jets 
spatially matched to partons (filled histogram). 



3. Parameterization of Momentum Resolution for Muons 

We describe the resolution of the central tracker 
through the uncertainty on the signed curvature of a 
track, the ratio of the electric charge and of the transverse 
momentum of a particle, parameterized as a function of 
pseudorapidity. The muon transfer function is parame- 
terized as 



Wft {k x , Ky) 



2ko- 



exp 



(16) 



where k x — (q/pT) x and K y = {q_/pT) y , with the charge 
q and transverse momentum px of the original muon (y) 
or its reconstructed track (x). The resolution 



2. Parameterization of Energy Resolution for Electrons 

The electron energy resolution is parameterized by the 
transfer function 



We (E x ,E y 



2iro- 



cxp 



1 fE x - EL 



where E x is the reconstructed electron energy, 
E' y = 1.000 • E y + 0.324 GeV, 
a= ^(0.028 ■ E' y f + {S ■ E' y f + (0.4 GeV) 
0.164 GeV 1/2 0.122 GeV 1/2 



(11) 

(12) 
(13) 



5 
C 



1.3519 



2.0956 GeV 6.9858 GeV 



E' 2 

v 



(15) 



E y is the energy of the original electron, and 9 e is the 
polar angle of the electron with respect to the proton 
beam direction. The parameters above are derived from 
the detailed modeling of electron energy response and 
resolution used in Ref. 13-411 . 



for Id < 1.4 
» , , . . ( 17 ) 



& 2 + {c-(\ri\-lA)y for \r]\ > 1.4 



is obtained from muon tracks in simulated events where 
the a and c parameters are linear functions of 1/pt- 



a = a +a 1 - 1/pr, 
c = c + ci • 1/pr- 



(18) 
(19) 



The values of the coefficients are given in Table IIIII 
for muon tracks with associated and no associated hits 
in the silicon tracker. This simplified parameterization 
of the momentum resolution is valid at high transverse 
momenta (pr > 20 GeV) where the limitations in coor- 
dinate resolution dominate over the effects of multiple 
scattering. 



VI. CALIBRATION OF THE MEASUREMENT 

The fully simulated MC samples described in Sec. II VI 
are used in ensemble studies to calibrate the result from 
the ME method by determining and correcting for biases 
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Table III. Parameters for muon transfer functions for muon 
tracks with and without hits in the SMT. 



Parameter 



With hits No hits 
in the SMT in the SMT 



cf (GeV- 1 ) 



1.125 x 1CT 



x 10" 



c (GeV- 1 ) 7.668 x 10" 3 2.070 x 10~ 



Cl 



1 

0.8 
0.6 
0.4 
0.2 




D0, e+jets (a) 



f = 0.76 x f 9 



0.2 0.4 0.6 0.8 1 

.gen 



1 

0.8 
0.6 
0.4 
0.2 




D0, u+jets (b) 



f = 0.81xf 9 



0.2 0.4 0.6 0.8 1 

.gen 



> 10 




■„p.05 



-0.05 



8 10 
172.5 GeV 




-0.05 



Figure 4. Mean values of (a) mt and (b) fejBS extracted from 
ensemble studies, as a function of the input values fitted to 
straight lines. Dashed lines represent 1:1 correlations of ex- 
tracted and input values. 



in the extracted parameters and their estimated uncer- 
tainties. Such biases can be due, for example, to limita- 
tions in the LO ME used in Eq. ((3]) or to the imperfect 
description of detector resolution using transfer functions 
with a limited number of parameters. Five tt MC samples 
are generated for mf n = 165, 170, 172.5, 175, and 180 
GeV, with two more produced from the 172.5 GeV sam- 
ple by re-scaling all jet energies by ±5%. P s i g and Pbkg 
are calculated for these samples and for the IV+jets MC 
samples. Events are drawn randomly from a tt sample 
with a particular mass and the IV+jets sample to form 
pseudoexperiments, each with a number of events equal 
to the one observed in data (before requiring > 1 6-tagged 
jets), with the signal fraction fluctuated according to a 
binomial distribution relative to that determined from 




Figure 5. Widths of the pull distributions for (a) m t and (b) 
fcjBS from ensemble studies as a function of the input values. 



Figure 6. Extracted signal fractions from ensemble studies 
as a function of the input values for the (a) e+jets and (b) 
/i+jets channels. 



data. The values of m t and &jes are extracted for each 
pseudoexperiment according to the procedure described 
in Sec. IV Al using only events with at least one 6-tagged 
jet. A thousand pseudoexperiments are performed for 
each of the 7 tt samples. The means (and their uncer- 
tainties) of all 1000 measured values of m t and fcjES m 
each sample are determined from Gaussian fits to their 
distributions and plotted versus the input mf on — 172.5 
GeV and fcjES _ lj respectively. A straight line is fitted 
to the plotted points, representing the response function 
used to correct the measurement from data (Fig. Q|. For 
each pseudoexperiment, we also calculate the pulls, de- 
fined as (m t ~ (m t )) / <r(m t ) and (fc JE s - (fcjEs))/o-(fcjEs), 
where (m t ) and (fcjEs) are the mean measured rrit and 
/cjes; respectively, for all pseudoexperiments, and cr(m t ) 
and cr(fcjEs) are the RMS of m t and fcjESi respectively, 
for the given pseudoexperiment. The width of the pull 
distributions for nit and fcjES are shown as a function 
of mf cn and fcj E g in Fig. [5l The average widths of the 
nit and fcj E s pull distributions are 1.08 and 1.07, respec- 
tively. 

The signal fraction for the ensemble studies is deter- 
mined from the selected data sample using the method 
described in Sec. IV Al To correct for biases in the de- 
termination of this fraction, a calibration is done using 
the IV+jets and 172.5 GeV tt MC samples, wherein 1000 
pseudoexperiments are performed using the same proce- 
dure as described in the previous paragraph, but with 
signal fractions set to a different value in each test. The 
extracted signal fractions as a function of their input val- 
ues are shown in Figs. Uta) and [U[b) for the e+jets and 
/i+jets channels, respectively. Straight lines are fitted 
to the points in plots representing the response functions 
used to correct the fractions determined from the selected 
data sample. The calibration of the signal fraction is per- 
formed separately for the e+jets and ^t+jets channels. 
The corrected fractions are 0.35 ± 0.05 and 0.41 ± 0.06 
for the e+jets and ^t+jets channels, respectively, prior to 
requiring at least one 6-tagged jet. These fractions are 
0.71±0.05 and 0.75±0.04 for the e+jets and ^+jets chan- 
nels, respectively, after requiring at least one 6-tagged jet. 
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VII. FLAVOR-DEPENDENT JET RESPONSE 
CORRECTION FOR MC EVENTS 

The validity of the calibration procedure described in 
the previous section is based on the assumption of a per- 
fect MC simulation of the events and of the detector. 
Uncertainties in this assumption are discussed in the sec- 
tion on systematic uncertainties (Sec. IIXj) . The in situ 
jet energy scale employed in this analysis can account 
for a global scale discrepancy between data and MC jet 
energy scales (see Sec. IIXB 1[) by rescaling the ener- 
gies of the two light jets from the W — > qq' decay of 
tt events to the world average mass of the W boson @ . 
This same rescaling is also applied to the two b jets in 
the event. However, jets originating from different par- 
tons have different kinematic characteristics and particle 
compositions. In particular, 6 and light jets with different 
electromagnetic fractions can lead to different responses 
in a non-compensating calorimeter. Such features, if not 
properly simulated, can result in a systematic shift in the 
determination of the top-quark mass. In fact, the largest 
contribution to the total systematic uncertainty of our 
previous analysis in Ref. || is the fe/light-quark response 
ratio which was an estimate of the effect of such a dis- 
crepancy. 

To bring the simulation of the calorimeter response 
to jets into agreement with data, and thereby reduce the 
systematic uncertainty associated with a jet response dif- 
ference in data and MC, we determine a flavor-dependent 
correction factor as follows. We note a discrepancy in the 
predicted energy deposition in the calorimeter between 
data and MC when we apply the single-particle responses 
from data and MC to the individual particles within MC 
jets that are spatially matched to reconstructed jets [35[ : 



Data 



V = 



J2 Ej ■ R 
Ei ■ Rf 10 ' 



(20) 



where the sums run over each particle i in the MC parti- 
cle jet, Ei is the true energy of particle i, and Rf ata and 
i?^ 10 are the single- particle responses in data and MC, 
respectively. We define a correction factor for a jet of fla- 
vor f3 (= light quark, gluon, or b quark) as the ratio of the 
discrepancy for jets of flavor j3 to the flavor-averaged dis- 
crepancy for jets in 7+jet events, F^ orl . = V 13 / (T> 1+ i ct y 
Defining the correction this way preserves the standard 
MC jet energy scale that is, strictly speaking, only ap- 
propriate for the 7+jet events from which it is derived. 
At the same time, it brings the relative response differ- 
ence between jets of flavor fj and jets in 7+jet events 
in MC into agreement with that in data. The quantity 
Ei OT r — 1 is shown in Fig. [7] as a function of jet pt and 
77 for light-quark, gluon, and b jets. The shaded band at 
F^„ — 1 = in each plot corresponds to the correction 
for jets in 7+jet events. We apply these correction fac- 
tors to the light-quark jets and b jets in a tt MC sample 
generated with mf en = 172.5 GeV, extract m t and /cjes 
using our analysis technique, and compare them with the 



values extracted from the same set of events without us- 
ing this correction. We find shifts of Am t = 1.26 GeV 
and AfcjEs = —0.005 relative to the uncorrected sample. 
Repeating this study on a tt MC sample appropriate for 
the previous analysis [§[ yields shifts of Am t = 1.28 GeV 
and A£:jes = -0.005. 



VIII. MEASUREMENT OF THE TOP-QUARK 

MASS 



The likelihoods L (x; m t ) and L (x; fcjEs) for the se- 
lected data, calculated according to Eq. © and Eq. (J7J, 
respectively, are calibrated by replacing m t and fcjES by 
parameters fitted to the response plots of Sec. I VI I 



= K- 172.5 ,GcV) - pT + m 5 GcVj 
Pi 

,,calib_ (fcjES-l)-Po' ES 
K-i 



MES 



fe.IES 

Pi 



+ 1, 



(22) 



where p™* and p^ JES are the parameters of the m t and 
fejES response functions shown in Fig.UJa) and Fig.fJJb), 
respectively, and mt and fcjES and their uncertainties are 
extracted from the mean and RMS values of the cali- 
brated likelihoods shown in Figs.[5ta) and(S]Jb). The ex- 
tracted uncertainties for m t and /cjes are multiplied by 
1.08 and 1.07, respectively, to correct for deviations of the 
average pull widths from unity (see Sec. IVI|) . Figure [S] 
shows the fitted Gaussian contours of equal probability 
for the two-dimensional likelihoods as a function of m t 
and fcjES- We find m t = 174.75 ± 1.28(stat + JES) GeV 
and fcjES = 1-018 ± 0.008(stat). Applying the shifts of 
Amt = 1-26 GeV and A/cjes = —0.005 described in Sec. 
IVIII yields a measured top-quark mass and jet energy 
scale factor of 

m t = 176.01 ± 1.28(stat + JES) GeV 

= 176.01 ± l.Ol(stat) ± 0.79(JES) GeV, 
fcjES = 1-013 ± 0.008 (stat). 

Distributions in expected uncertainties, determined from 
1000 pseudoexperiments performed on the MC tt sample 
for mf en = 175 GeV, are shown in Figs. fTCT a) and [TOT b) 
for mt and ZcjeSi respectively. The measured uncertain- 
ties, indicated by the arrows, are within the expected 
range observed in MC, and do not depend in any appre- 
ciable way on the assumed value of m t . 



IX. SYSTEMATIC UNCERTAINTIES 

We evaluate systematic uncertainties for three cate- 
gories. The first category, modeling of production, ad- 
dresses uncertainties in the MC modeling of tt and 
W+jets production. The second category, modeling of 
detector, deals with the uncertainties in jet energy and 
lepton momentum scales and the simulation of detector 
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Figure 7. (color online) Correction factors for data-MC jet response difference for light-quark, gluon, and 6 jets as a function of jet pt and r\. Light shaded bands 
represent statistical uncertainties and dark shaded bands represent systematic uncertainties 



14 



E 1 



D0, 2.6 fb 1 (a) 
m t = 176.01 +1.28 GeV 
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m, GeV 

Figure 8. Calibrated projections of the data likelihoods onto 
the (a) mt and (b) A:jes axes with 68% confidence level re- 
gions indicated by the shaded areas. The values of mt and 
fcjES shown in the figures are after applying all the corrections 
described in Sec. IVIII and Sec. I Villi 
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response and any associated efficiencies. The third cate- 
gory involves uncertainties in the calibration of both m t 
and the signal fraction /, and a possible bias from the 
exclusion of multijet events in MC ensemble studies. The 
contributions to the systematic uncertainty are summa- 
rized in Table IIVI In the first three sections below, we 
describe the evaluation of each of these contributions in 
more detail. In the fourth section, we discuss how sys- 
tematic uncertainties from the previous analysis of 1 fb -1 
of integrated luminosity Q are updated to facilitate the 
combination of the two results presented in Sec. IXII Ex- 
cept for the Data-MC jet response difference described in 
Sec. IIX B 21 all of the systematic uncertainties described 
below are calculated prior to the flavor-dependent jet re- 
sponse corrections of Sec. IVIII 

We adopt the following convention for systematic un- 
certainties 6m t in m t , and classify them into two types. 
The first type, referred to as the Type I uncertainty, is the 
effect of the ±1 sd variation of a relevant quantity. The 
second type, referred to as the Type II uncertainty, is due 
to the difference between models. For Type I uncertain- 
ties, we refer to the central or default value of the mea- 
surement as m° and to the measurement corresponding 
to the 1 sd (—1 sd) variation as mf (m^~). We compute 

\ m t - m t\l 2 
and according 



Type I uncertainties according to Smt 



m f < 



mt 
\ m t 



or mf 



< mt < m 



A 



m 



/2 



if 



have the same sign. We compute Type II 



if m t < 
to 8m t = 
and 

uncertainties by taking the maximal difference between 
the models as the + and — systematic variations. 

Many of our systematic uncertainties are evaluated by 
comparing two MC tt samples generated with the same 
input mass mf cn . For these studies, we use samples with 



a value of : 



,gon 



close to the world average of m t such as 



172.5 GeV or 170 GeV. 



A. Modeling of Production 



Figure 9. (color online) Fitted contours of equal probability 
for the two-dimensional likelihood L (x; mt, fcjEs)- 




1.5 
a(rn) GeV 



Figure 10. Expected uncertainty distributions for (a) mt and 
(b) fcjES determined from 1000 pseudoexperiments performed 
on the MC tt sample for mf cn = 175 GeV. The measured 
uncertainties are indicated by the arrows. 



1. Higher- Order Effects 

The MC tt samples used to calibrate our measurement 
are generated using ALPGEN for the hard-scattering pro- 
cess and PYTHIA for shower evolution and hadronization 
(Sec. HV[) . We compare the LO generator ALPGEN with 
the next-to-leading order MC generator MC@NLO [3(| . 
in order to evaluate possible contributions from higher- 
order effects such as additional radiation of hard jets or 
gg contributions. We compare ALPGEN and MC@NLO 
MC tt samples with identical values of mf cn that both 
use HERWIG for shower evolution and hadronization. 
HERWIG is used in both cases for consistency because 
MC@NLO can only be used with HERWIG (alpgen can be 
used with PYTHIA or HERWIG) and we are not interested 
in comparing different models for shower evolution and 
hadronization in this study. Ensemble studies are per- 
formed on both samples and the difference in the mean 
extracted m t from ensembles for the two samples is found 



15 



Table IV. Summary of systematic uncertainties. 



Source 



Uncertainty (GeV) 



Modeling of production: 
Modeling of si 



Higher-order effects 


±0. 


2"> 


ISR/FSR 


±0, 


26 


Hadronization and UE 


±0 


58 


Color reconnection 


±0 


28 


Multiple pp interactions 


±0 


07 


Modeling of background 


±0. 


16 


W +jets heavy- flavor scale factor 


±0 


07 


Modeling of b jets 


±0. 


09 


Choice of PDF 


±0 


24 


Modeling of detector: 






Residual jet energy scale 


±0. 


21 


Data-MC jet response difference 


±0. 


28 


h-T a crcri n cr pffipipnpv 


±0, 


08 


Trigger efficiency 


±0, 


01 


Lepton momentum scale 


±0, 


17 


Jet energy resolution 


±0, 


32 


Jet ID efficiency 


±0, 


26 


Method: 






Multijet contamination 


±0, 


14 


Signal fraction 


±0, 


10 


MC calibration 


±0, 


20 


Total 


±1. 


02 



3. Hadronization and Underlying Event 

In simulating partem evolution and hadronization, 
PYTHIA and HERWIG model the parton showering, 
hadronization, and underlying event (UE) differently. 
To estimate the impact of this difference, we compare 
two MC tt samples with identical values of mf en , using 
ALPGEN for the hard-scattering process, but one sam- 
ple using pythia and the other using HERWIG for par- 
ton showering and hadronizaton. Ensemble studies indi- 
cate a difference in the means of the extracted mt to be 



4- Color Reconnection 

The MC samples used in this analysis do not simulate 
color reconnection for the final-state particles [3{j. To 
evaluate the possible effect of color reconnection on the 
determination of m t , we compare two MC tt samples with 
identical values of mf on , using PYTHIA 6.4 tunes Apro 
and ACRpro, which are identical except for the inclusion 
of color reconnection in ACRpro. Ensemble studies of tt 
events performed on both samples yield a difference in the 
means of the extracted m t of m* pro — m t CRpro = 0.26 ± 
0.28 GeV. We take the uncertainty on this difference and 
assign ±0.28 GeV as the contribution from this source. 



5. Modeling of Jet Mass 



to be m t MC @ NL ° _ m£ LPGEN = 0.10 ± 0.25 GeV. Here, as in 
all the other systematic sources described below, when a 
shift in the value of the estimated parameter is statisti- 
cally dominated, we replace the shift with its statistical 
uncertainty for the estimate of uncertainty. We, there- 
fore, assign an uncertainty of ±0.25 GeV as the contri- 
bution from this source. 



2. ISR/FSR 

The uncertainties from this source are in the modeling 
of additional jets due to initial and final-state radiation 
(ISR/FSR). To evaluate this contribution, we compare 
three pythia samples having identical values of mf en , 
with input parameters taken from a CDF ISR/FSR study 
based on the Drell-Yan process |38j . The three sets of 
parameters correspond to a fit to data and ±1 sd excur- 
sions. Half of the difference between the two excursions 
corresponds to a change in m t of 0.26 ± 0.19 GeV. 



Unlike the jet algorithm used in Run I of Tevatron, the 
iterative midpoint cone algorithm used for Run II defines 
jets of intrinsic mass [2l| . The effect of inaccuracies in the 
simulation of jet masses on the top-quark mass measure- 
ment is found to be negligible and is presently ignored. 



6. Multiple pp Interactions 

Effects from additional pp interactions are simulated by 
overlaying on MC events unbiased triggers from random 
pp crossings. These overlaid events are then reweighted 
according to the number of interaction vertices to assure 
that the simulation reflects the instantaneous luminosity 
profile of the data. To evaluate the contribution from the 
uncertainty associated with the reweighting procedure, 
we repeat the ensemble studies used to derive the mt 
calibration, but without the reweighting. The rederived 
calibration is applied to L (x; mt) for the selected data 
sample, m t is extracted and compared with the value 
from the default calibration, and found to shift by —0.07 
GeV. This extreme check of the size of this contribution 
to the uncertainty shows that our result is not affected 
significantly by variations in luminosity. 
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7. Modeling of Background 

This systematic uncertainty receives contributions 
from two sources, one based on the data-MC discrep- 
ancy in background-dominated distributions, and a sec- 
ond from uncertainty in the renormalization scale used 
to generate the JV+jets samples. For the first source, we 
identify distributions in which there is poor agreement 
between data and MC in the modeling of background. 
Specifically, in both channels, we examine lepton px and 
the n of the jet of lowest px in the 3-jet multiplicity bin. 
Ensemble studies are performed on a sample of MC tt 
events using background events reweighted to match the 
distributions in data. The mean of the extracted nit for 
this sample is found to shift by —0.03 GeV relative to that 
of the same MC tt events using the default background 
events. 

The TV+jets MC samples used in this analysis (Sec. 
IIVP are generated using identical renormalization and 
factorization scales of fj, = Myy + ^ p T where the sum is 
over the jets in an event. To evaluate the effect of the un- 
certainty in this scale, we generate two more JV+jets MC 
samples with modified renormalization and factorization 
scales of /i/2 and 2/z. We perform ensemble studies on a 
tt MC sample using these modified IV+jets samples, and 
find that the means of the extracted m t shift by 0.13 GeV 
(///2) and 0.32 GeV (2fi) relative to the studies using the 
default tV+jets sample. We take half of the larger excur- 
sion and assign ±0.16 GeV as the contribution from this 
source. 

The contributions from the above data-MC discrep- 
ancy for the background and from the uncertainty on 
the scales are combined in quadrature for a total of sys- 
tematic uncertainty of ±0.16 GeV. 

8. W+jets Heavy-Flavor Scale Factor 

The default heavy-flavor content in LO ALPGEN MC 
TV+jets (Sec. IIV|) is increased by a factor of 1.47 for the 
M / cc±jets and Wbb+jets contributions to achieve agree- 
ment with NLO calculations of cross sections that include 
NLL corrections based on the MCFM MC generator (2(|. 
To evaluate the uncertainty from this source, we shift 
this factor up to 1.97 and down to 0.97 and, for each 
variation, repeat the ensemble studies described in Sec. 
IVII for the calibration of nit, apply this to L (x; m t ) in 
data, and re-extract nit- The shifts in m t relative to the 
default value are found to be —0.07 GeV and 0.02 GeV 
when the scale factors are shifted up and down, respec- 
tively. We assign ±0.07 GeV as the contribution from 
this source to the uncertainty of nit ■ 

9. Modeling of b jets 

Possible effects in modeling 6-quark fragmentation are 
studied by reweighting the simulated tt events used in the 



calibration of the measurement to simulate other choices 
of &-quark fragmentation models for the b jets. All the 
default MC samples used in this analysis consist of events 
that are reweighted from the default pythia 6-quark 
fragmentation function (based on the Bowler model (io| ) 
to a Bowler scheme with parameters tuned to data col- 
lected at the LEP e + e~ collider [4l|. To evaluate the 
systematic uncertainty, these events are reweighted again 
to account for differences between LEP and SLAC e + e~ 
data [4l[ . The ensemble studies of m t are repeated using 
these reweighted events, the new calibration applied to 
L (x; m t ) for data, and m t extracted, nit is found to shift 
by 0.08 GeV relative to the default value. 

Additional differences in the response of b jets can be 
expected in the presence of semileptonic decays of b or 
c-quarks. The incorrect simulation of semileptonic b and 
c-quark decay branching fractions can therefore lead to 
a systematic shift in the extracted value of nit- We take 
an uncertainty of ±0.05 GeV determined in Ref. [42| as 
the contribution from this source. 

Combining the two above uncertainties in quadrature 
gives ±0.09 GeV, which we assign as the systematic un- 
certainty for the modeling of b jets. 



10. Choice of PDF 

We evaluate this systematic uncertainty using a 
pythia MC tt sample that is reweighted to match possi- 
ble excursions in the PDF parameters represented by the 
20 CTEQ6M uncertainty PDFs [13 . Ensemble studies 
are repeated for each of these variants for only tt events, 
and the uncertainty evaluated using the following formula 

1 / 20 \ 1 / 2 

Smf DF = g \Y[AM(S+) - AM(S~)} 2 ) (23) 

where the sum runs over PDF excursions in the positive 
(5+) and negative (S~) directions. 5mf DF is found to 
be 0.24 GeV. 



B. Modeling of Detector 

1. Residual JES Uncertainty 

The in situ jet energy calibration employed in this anal- 
ysis addresses a possible global scale difference in JES 
between data and MC. Any other discrepancy, such as a 
dependence on px and n, can have a systematic effect on 
the determination of nit- To estimate this, the fractional 
uncertainty associated with the standard jet energy cor- 
rection, derived using the 7±jet and dijet samples, is 
parameterized as a function of px and n. This uncer- 
tainty includes statistical and systematic contributions 
from both data and MC added in quadrature. All jet 
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energies in a tt MC sample are then scaled up by the pa- 
rameterized uncertainty as a function of and r\. The 
parameters are then shifted in such a way that the av- 
erage scale shift applied to all jets vanishes. Ensemble 
studies are performed on the default and scaled samples, 
and the extracted mt found to shift by 0.21 GeV relative 
to the default sample. 



2. Data-MC Jet Response Difference 

The uncertainties in the flavor-dependent jet response 
correction for MC events (described in Sec. IVII|) . used to 
bring the simulation of calorimeter response into agree- 
ment with that observed in the data, are associated with 
uncertainties in single-particle responses in data and MC. 
To evaluate the effect of these uncertainties on the value 
of m t , we change the correction factors by ± 1 sd and ap- 
ply them to the light jets and b jets in a ti MC sample. 
The value of mt is extracted and the mean is found to 
shift by ±0.28 GeV relative to the sample corrected using 
the central values. 



3. b-Tagging Efficiency 

Discrepancies in the 6-tagging efficiency between data 
and MC can lead to a systematic shift in the extracted 
m t . To evaluate the effect of possible discrepancies, the 
tag rate functions for b and c quarks and the mistag rate 
function for light quarks are changed by 5% [23[ and 20%, 
respectively, corresponding to the uncertainties on these 
functions. Ensemble studies for all tt MC samples are 
then repeated and the mt calibration rederived and ap- 
plied to data to extract m t . The result is compared with 
that from the default calibration and found to shift by 
-0.08 GeV. 



4- Trigger Efficiency 

The MC events used in this analysis have associated 
weights to simulate the effect of trigger efficiencies. To 
evaluate the effect of the uncertainties in these weights on 
the top-quark mass, we repeat the ensemble studies on 
all tt MC samples with the weights set to unity, rederive 
the mt calibration, and apply it to the data to extract 
m t . The result is found to shift by —0.01 GeV. 



5. Lepton Momentum Scale 

A relative difference in the lepton momentum scale be- 
tween data and MC can have a systematic effect on mt- 
To evaluate this, we first determine the size of the dis- 
crepancy and correct the scale of one ti MC sample. En- 
semble studies are repeated on the corrected sample and 



the mean of the extracted mt is found to shift by 0.17 
GeV relative to the default sample. 

6. Jet Energy Resolution 

Since the jet transfer functions used are derived from 
MC samples, improper simulation of jet energy resolution 
can result in a bias in the extracted m t . To evaluate a 
possible bias, ensemble studies are performed using a tt 
MC sample with jet energy resolutions degraded by 1 sd. 
The mean of the extracted mt in this sample is found to 
shift by 0.32 GeV Q. 

7. Jet ID Efficiency 

The uncertainties associated with the scale factors used 
to achieve data-MC agreement in jet ID efficiencies are 
propagated to the measurement of m t by decreasing the 
jet ID efficiencies in a ti MC sample according to these 
uncertainties. We can only simulate a decrease and not 
an increase, as reconstructed jets can be dropped but not 
created. Ensemble studies indicate that the mean of the 
extracted m t shifts by 0.26 GeV relative to that of the 
default sample |44ll45j. 

C. Method 

1. Multijet Contamination 

The multijet background is not included in the ensem- 
ble studies used to derive the calibrations described in 
Sec. I VII as we have assumed that Pbkg 3> -Fsig for such 
events (see Sec. IV A[) , resulting in a negligible influence 
on the determination of m t . To evaluate possible system- 
atic effects due to this assumption, we select a multijet- 
enriched sample of events from data by inverting the lep- 
ton isolation criterion in the event selections. We repeat 
the ensemble studies to derive the mt calibration using 
the multijet-enriched sample in the sample composition. 
The rederived calibration is applied to data and the ex- 
tracted m t is found to shift by 0.14 GeV relative to the 
default calibration (45| . 

2, Signal Fraction 

The signal fractions determined from data and used 
in the ensemble studies have associated statistical uncer- 
tainties. These signal fractions are varied by their uncer- 
tainties, independently for each decay channel, and the 
ensemble studies repeated for all MC samples to rederive 
the mt calibration shown in Fig. |4ja). The new calibra- 
tions are then applied to the data and results compared 
with those obtained using the default calibration. The 
resulting uncertainties in m t evaluated by changing the 
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signal fractions in each decay channel are then added in 
quadrature and divided by two to obtain a total of ±0.10 
GeV. 



3. MC Calibration 

We estimate the effect of the statistical uncertainties 
associated with the offset and slope parameters deter- 
mined from the fit to the response plot shown in Fig. 
|4ja). To estimate this uncertainty, we change these two 
parameters, one at a time, by their uncertainties, and 
apply the modified calibration to the data to extract m ( , 
and calculate the difference relative to the mt extracted 
using the default calibration. We combine, in quadra- 
ture, the differences in m t resulting from such changes in 
each parameter, and find an uncertainty of ±0.20 GeV. 



D. Treatment of Systematic Uncertainties in 
Previous Analysis 

To facilitate the combination of the new measurement 
with the previous one, we have updated the system- 
atic uncertainties presented in Table I of Ref. All 
of the uncertainties in this table are unchanged, except 
for the uncertainties in the modeling of signal and the 
relative 6/light-quark response ratio. The uncertainty 
for the modeling of signal in the previous analysis is 
replaced with one from the current analysis, which in- 
cludes contributions from uncertainties in the modeling 
of higher-order effects, ISR/FSR, hadronization and un- 
derlying event, color reconnection, and multiple hadron 
interactions. The uncertainty on &/light-quark response 
is replaced with that associated with differences in jet 
response in data and MC for the current analysis (see 
also Sec. IVII[) . The uncertainty in the modeling of back- 
ground in Table I of Ref. Q is the sum in quadrature of (i) 
the uncertainty in the heavy-flavor scale factor, and (ii) 
the uncertainty associated with discrepancies between 
data and MC background distributions. Since the un- 
certainty on the renormalization and factorization scale 
was not evaluated in the previous analysis, we include 
the additional contribution described in the second part 
of Sec. IIX A 71 We also evaluate the uncertainty asso- 
ciated with the flavor-dependent jet-response correction 
factors appropriate for the previous analysis, using the 
procedure described in Sec. IIX B 21 We find the mean of 
the extracted m t shifts by 0.13 GeV (—0.22 GeV) relative 
to the sample corrected with the central values when we 
change the correction factors by 1 sd (— 1 sd). We assign 
±0.22 GeV as the contribution from this source. Adding 
the contributions from all sources in quadrature gives a 
total of ±0.97 GeV. 



X. RESULT OF THE CURRENT 
MEASUREMENT 

We measure the mass of the top quark in tt lepton+jets 
events using a matrix element method that combines an 
in situ jet energy calibration with additional information 
from the standard jet energy scale derived from 7+jet 
and dijets samples. Using data corresponding to 2.6 fb _1 
of integrated luminosity collected by the DO experiment 
from Run II of the Tevatron collider, we extract the value: 

m t = 176.01 ± l.Ol(stat) ± 0.79(JES) ± 1.02(syst) GeV, 

or m t = 176.01 ± 1.64 GeV. 



XI. COMBINATION WITH THE PREVIOUS 
MEASUREMENT 

Our result from a previous measurement using the 
same analysis technique, and based on earlier data cor- 
responding to 1 fb _1 of integrated luminosity, is mt = 
171.5 ± 1.76(stat + JES) ± l.l(syst) GeV Q. Applying 
the shift of Am t = 1.28 GeV described in Sec. 17111 
and using updated systematic uncertainties described in 
Sec. ITXDl yields 

mt = 172.74 ± 1.44(stat) ± 1.05(JES) ± 0.97(syst) GeV, 

or m t = 172.74 ±2.03 GeV. 

We combine the two measurements using the BLUE 
method (4(| H3| to get a result equivalent to 3.6 fb _1 of 
integrated luminosity. The combined value of the mass 
is 

mt = 174.94 ± 0.83(stat) ± 0.78(JES) ± 0.96(syst) GeV, 

or m t = 174. 94± 1.49 GeV. The procedure we follow uses 
the same method and classes of uncertainty as used by 
the Tevatron Electroweak Working Group [5j in combin- 
ing individual measurements for Tevatron averages of the 
top-quark mass. 
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